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Abstract 

We consider the early universe at temperatures close to the fundamental scale of 
gravity {Md Mpianck) in models with extra dimensions. At such temperatures a 
small fraction of particles will experience transplanckian collisions that may result in 
microscopic black holes (BHs). BHs colder than the environment will gain mass, and 
as they grow their temperature drops further. We study the dynamics of a system (a 
hlack hole gas) defined by radiation at a given temperature coupled to a distribution of 
BHs of different mass. Our analysis includes the production of BHs in photon-photon 
collisions, BH evaporation, the absorption of radiation, collisions of two BHs to give a 
larger one, and the effects of the expansion. We show that the system may follow two 
different generic paths depending on the initial temperature of the plasma. 



1 Introduction 



Large (ADD) [IJ or warped (RS) [2J extra dimensions open the possibility that the funda- 
mental scale of gravity Mo is much lower than Mp ^ 10^^ GeV. This could imply that the 
transplanckian regime [3] is at accesible energies. Collisions in that regime are very different 
from what we have experienced so far in particle colliders: due to its spin 2 the graviton 
becomes strongly coupled and dominates at distances that increase with the center of mass 
energy y/s. In particular, at small impact parameters one expects that gravity bounds the 
system and the two particles collapse into a microscopic black hole (BH) of mass around 
y^. Such BH would evaporate [4] very fast (the time scale is set by 1/Md) into final states 
of high multiplicity, a possible LHC signature that has been extensively discussed in the 
literature [Sj. 

These collisions, however, may have also occurred in the early universe if the temperature 
was ever close to M^- Consider a period of inflation produced by a field on the brane, 
followed by reheating into 4-dimensional (4-dim) species. If the reheating temperature is 
Trh > O.lMp), particles in the tail of the Boltzmann distribution may collide with enough 
energy to form BHs. Hot BHs will evaporate, but if a mini-BH is colder than the environment 
it absorbs more than it emits growing and becoming colder, which in turn increases its 
absorption rate (see below). The fact that heavier BHs live longer distinguishes them from 
massive particles or string excitations possibly produced at these temperatures, since the 
lifetime of the latter is inversely proportional to their mass. As a consequence, one expects 
that BHs are a critical ingredient at temperatures near the fundamental scale. Notice that 
these BHs are not the primordial ones formed from the gravitational collapse of density 
fluctuations [71 El El HO] at lower temperature T. 

Although BHs would also be present in a 4-dim universe at T ~ Mp, they are most 
relevant in TeV gravity models. The reason is that the expansion of the universe is a long- 
distance process, so its rate is dictated by Mp (not Mp,) also in models of TeV gravity. If the 
bulk is basically empty [HI [121 IE] or thin (see below), then the expansion at T < Mp, will 
be negligible in terms of the fundamental time scale 1/Mp, and there is plenty of time for 
collisions producing BHs to occur and for the BHs to grow. In contrast, in a 4dim universe 
the expansion rate at T ^ Mp is of order H^^ ^ 1/Mp, the temperature of the universe 
drops before BHs have grown, and once it goes below TpH they evaporate in the same time 
scale. 

In this article we explore the implications of having initial temperatures near the funda- 
mental scale of gravity. First we define a consistent set up for TeV gravity. Then we study 
the behaviour of a single BH inside a thermal bath in an expanding universe. Finally we 
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consider the generic case, a black hole gas, with radiation coupled to a distribution of BHs 
of different mass. We find remarkable that the effect of these mini-BHs in the early universe 
has been almost completely overlooked in the literature (the only analysis that we have 
found is given by Conley and Wizansky in [M]), although many authors have considered 
temperatures close (and above) the Planck scale (see [15] and references therein). 

2 Consistent TeV gravity models 

There are two generic frameworks that may imply unsuppressed gravity at ^ 1 TeV (or 
at any scale M^) < Mp). In the first one (ADD), n compact dimensions of length L < 1 mm 
introduce a large number of KK excitations of the graviton. These gravitons, of effective 
4-dim mass proportional to rric = L~^, couple very weakly ^/s/Mp) to ordinary matter. 
However, the large number of effective gravitons ~ {\fs/mcY involved in a collision gives an 
amplitude of order one. 



In contrast, in the second scenario (RS) the KK excitations of the graviton have unsuppressed 
couplings to matter but large masses, right below Mc, so a few KK modes suffice to define 
an order one gravitational interaction at that scale. 

The ADD set up has the basic cosmological problem pointed out in [12]. Essentially, at 
temperatures T ^ ttIc KK gravitons will be abundantly produced in annihilations of brane 
particles, due to the large multiplicity of final states. If the initial temperature is large these 
gravitons will change the expansion rate at the time of primordial nucleosynthesis. Even if 
the initial temperature is as low as 1 MeV, their late decay will distort the diffuse gamma ray 
background in an unacceptable way. Obviously, an initial temperature close to Mp, would 
bring too many massive gravitons. 

One solution would be to consider RS models of TeV gravity. There, at T Mp, 
bulk and brane species have similar abundances, but massive gravitons will decay fast once 
T < nic ~ O.IM/), returning all the energy to the brane and defining acceptable 4-dim 
cosmological models. 

Even within the ADD framework, however, we can consider hybrid models where the 
connection between ttIc and M^) is not the one given in Eq. ([2]). This can be obtained, for 




(1) 



at ^ = Md if 
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example, with a warp factor [16]. The effect would be to push the KK modes towards the 
4dim brane, reducing the effective compact volume to V ~ (l/mc)" while increasing their 
coupling to matter, 

^j_(^y . (3) 



Ml M2 \Md, 

In this way, a smaller number of KK modes will imply an order one gravitational interaction 
at the same scale a/s = M^. If the free parameter rric takes the value in Eq. ([2]) we recover 
ADD, whereas for rric approaching Md we obtain RS. At distances below l/rric gravity would 
be higher dimensional (similar to ADD) whereas at larger distances it becomes 4-dimensional 
(like in the usual RS scenario). 

In this framework the KK modes of the graviton are not produced at temperatures 
T < rric. Therefore, > 10 MeV avoids astrophysical bounds [17] for any number n of 
extra dimensions [16] and Md = 1 TeV. On the other hand, as these KK gravitons have 
stronger couplings to matter, they can decay much faster than in the usual ADD model 
or decouple at temperatures below their mass. By changing rric (with Md fixed) it seems 
easy to obtain models with no gravitons at the time of primordial nucleosynthesis that are 
consistent with all cosmological observations. 



3 Single black hole in a thermal bath 

Let us consider the hybrid framework described in the previous section with n extra dimen- 
sions, a fundamental scale Md, and an independent compactification mass rric > 1 GeV. 
Gravity at energies above Md and distances smaller than l/rric is strongly coupled and 
higher dimensional (just like in ADD), so the radius and temperature of a BH of mass M 
are 

an /My/("+i) ^ n+1 

with ^ 

• (5) 

n+2 J ^ ' 

Once a BH reaches a radius ru = ^ I rric and fills up the whole compact space its (4-dim) 
size will not keep growing, since all the KK gravitons but the zero mode provide short 
distance interactions. The radius should start growing significantly only when the usual 
4-dim horizon (produced by the massless graviton) is of order 1/mc, i.e., for BH masses 
M ^ Mp/rric. Above this threshold BHs are basically 4-dim. 
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We will assume that the particles in the brane and the bulk (with g^, and degrees of 
freedom, respectively) are initially in thermal equilibrium at a temperature T = Tq. Notice 
that reheating in just the brane after a period of inflation would not justify an empty bulk, as 
in a time of order M'jy/T^ most of the energy would escape into KK modes. Here, however, 
the bulk may be much thinner and emptier than in the usual ADD model, implying a slow 
expansion in terms of the fundamental time (see below). 

It is easy to see that if the plasma temperature T and I/th are larger than ttIc a BH at 
rest in the brane will change its mass according to [H] 

" ^""^^^'^g. [t' - + g,c^ i^T^- nn)) , (6) 



where = g^TT'^/120, c„ is an order 1 coefficient that depends on the geometry of the 
compact space, and we have neglected gray-body factors [18]. The expression above assumes 
a thermalized photon (in the brane) and graviton (in the bulk) plasma, so it requires that 
changes occur slowly. In particular, if the BH grows very fast one has to make sure that the 
gain in mass is always smaller than the total energy of the plasma in causal contact with 
the hole: 

dM 4^ / T"\ 
['*^'''-^) ' 

where we have assumed t > l/m^. 

As the BH grows it enters a new phase when its mass reaches 

M,^Mn(^y^\ (8) 

which corresponds to a radius rn = ^/ttLc filling up the whole extra volume. For M > Mi 
the BH keeps gaining mass as far as the plasma temperature T is above Tbh ~ itt-c- However, 
as explained above, KK fields do not reach distances beyond vh ~ ^/itIc so the BH radius 
(and its temperature) will stay basically constant. 

Finally, masses above M2 ~ Mp/rric turn the BH into a purely 4 dimensional object: its 
radius 

2M 

grows larger than l/vc and the BH becomes too cold to emit massive gravitons. In this 
regime, if T is larger than rric the BH changes its mass according to 

dM vr / yn\ 2-4 
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whereas at lower plasma temperature it goes as 

dM TT / 

BH 



Eq.([H]) implies that lighter (hotter) BHs evaporate with an approximate lifetime 



1 / M 



3+n 
1 + n 



(12) 



and that BHs colder than the plasma will gain mass: heat flows from the hot plasma to the 
cold BH, but the effect is to cool the BH further and increase T — Tbh- This is, indeed, a 
very peculiar two-component thermodynamical system. 

On the other hand, the expansion of the universe is also affected by the presence of 
bulk species. We will assume that the extra dimensions are frozen (do not expand) and 
will integrate the matter content in the bulk. The large values of that we will consider 
imply large enough couplings with brane photons, so that the (equilibrium) abundance of 
KK modes of mass larger than the plasma temperature will be negligible. Therefore, at 
T > nic we have a universe with a radiation density 

Prad^^T \9*+gbCn—\ , (13) 

and an expansion rate 

If SttG 

= -^Prad , (14) 

i.e., Prad = Prado{Rol RY- Noticc that the second term in Eq. (fT3|) will slow down the change 
in the plasma temperature T{t) due to the expansion. In particular, if the bulk energy 
dominates then T oc ^2/(4+") times larger than a Hubble time. At temperatures below nic 
this term vanishes exponentially and all the bulk energy is transferred to the brane. 

To ilustrate the orders of magnitude involved, let us discuss a toy model with M/) = 1 
TeV, n = 1 (ci = 0.3), and just photons ((?* = 2) and gravitons {gb = 5) at Tq = 100 GeV. 
A BH of mass M < Merit = 12 TeV would be hotter than the environment, and it would 
evaporate in a time of order r ^ 10^'^ GeV~^ = 6.5 x 10"^*^ s. The Hubble time of a universe 
at this temperature is 

H-^ = Z = 1.4 X 10^^ GeV"^ = 9.2 x 10"^ s . (15) 
R 

A BH of initial mass Mq = 100 TeV will have a starting temperature of 8.7 GeV. In a time 
of order 17 GeV~^ its mass is already around 4.7 x 10'' GeV and its radius as large as the size 
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Figure 1: Evolution of a single BH of mass Mq = 100 TeV for Md = 1 TeV and n = 1 in an 
expanding universe at To = 100 (lower line) and 200 GeV (upper line). 

of the extra dimension (0.1 GeV^^). Then the BH keeps growing at approximately constant 
rate, since its size and temperature {Tbh ~ 1-6 GeV) change little with the mass. In a 
Hubble time the BH reaches a mass M ^ 5.7 X 10^*^ GeV. At later times the expansion cools 
the plasma, which slows the growth of the BH. Its maximum mass, M ^ 2 x 10^^ GeV, is 
achieved when T = Tbh at times of order t ^ 10^^ GeV~^. Finally, the BH will evaporate 
after r ^ 10^^ GeV^^ ^ 1 s. 

In Fig. 1 we plot this case together with the mass evolution for a larger initial temperature. 
To = 200 GeV. In this second case the expansion rate is faster (the Hubble time ^ 
2.6 X 10^^ GeV^^ is shorter), but the higher radiation density provides larger BH masses. 
Lower (higher) values of (M/j) would also imply larger BH masses. 

In none of the two cases described above the BH becomes purely 4 dimensional, i.e., with 
a mass M > 10^^ GeV. These large values of M can be obtained increasing the number n of 
extra dimensions and/or the ratio Mu/rric. A BH radius larger than l/m^, however, would 
not stop its growth, on the contrary, it would enhance the absorption rate of radiation by 
the BH. What stops the growth of a BH in any TeV-gravity framework is just the drop in 
the temperature of the plasma due to the expansion. 

It is also important to emphasize that the qualitative features of the process described 
above do not depend much on the details of the compactification (the length or the shape) 
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of the extra dimensions. They just depend on the fact that the fundamental scale of gravity 
Md is low and that the initial temperature of the radiation is close to it. In particular, the 
sequence of events would be similar in RS models, although there is close to M^) and the 
BHs are smaller [th ~ l/'^^c, and the absortion rate grows with the BH area). 

Finally, notice that the usual cosmological constraints [T3] on the initial (reheating) 
temperature of the early universe or the astrophysical bounds on are actually a probe 
of 1/L = mc, so they are avoided if irtc is above 1 GeV. 



4 Dynamics of a black hole gas 

Let us now deduce the equations that describe the production and evolution of BHs inside 
a thermal bath of temperature T. If T < Mc, < M the BHs will be produced in collisions 
between radiation in the high-energy tail of the Boltzmann distribution /^(A;). Once formed 
these BHs will be non-relativistic, with a kinetic energy K ^ T negligible versus the mass M 



and a velocity v ~ y 2T/M. We will then assume that the two components of the system are 
well described by the temperature T(t) of the plasma and the distribution /(M, t) expressing 
the number of BHs of mass M per unit mass and volume at time t. The total energy density 
is 

71-2 / T"\ /■ 

P{t) = Pradit) + PBH{t) ~ ^ [d* + 9bCn—\ + J dM M f{M,t) . (16) 



We identify the following processes changing the number density of BHs of mass M. 

• BH production in photon-photon collisions. Photons will collide with a cross section 
cr ^ irrfj to form a BH of mass M ^ y/s: 

(^^^) H^l^ Uh)Uk,) a{M) \v,-v,\ 5(vWT^-M) , 

(17) 

with Vi = 1. If T ^ M this expression can be approximated in terms of a modified 
Bessel function, 



To simplify our analysis we will neglect BH production in collisions of bulk particles. 
Although KK modes dominate the energy density, their cross section to form a BH is 
smaller (their wave function is diluted within the bulk), so this would be an order one 
contribution. 
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The collision of two BHs, of mass Mi and M2, to form a BH of mass M = Mi + M2: 



( dfiM,ty 

[ dt , 



dMidMs f{Mi,t)f{M2,t) (t{Mi,M2) vi2 6{Mi + M2-M) , 

(19) 



where the BH velocity is Vi = yST/Mj and vu = {\vi — V2\). We will take a BH-BH 
cross section of 

a{Mi,M2) = 7r{rHi+rH2f . (20) 

Notice also that for a minimum BH mass of Md, this contribution to /(M, t) is nonzero 
only if M > 2Md. 

A BH of mass M may collide with any other BH, which would reduce f{M,t): 

'dMi f{M,t)f{Mi,t) a{M,Mi) voi , (21) 



( dfjM^ty 
[ dt 



MMi->M2 



We can describe the absorption and emision of radiation using dM/dt in ([6]). The BHs 
of mass M will have a mass M + dM at t + dt, i.e., 

f{M,t) = f{M + dM,t + dt) . (22) 

This implies 

( df{M,t) \ _ df{M,t) dM 

[ dt L~ dM dt ■ 

\ / abs/em 

Finally, we have to add the effect of the expansion. The 4-dim scale factor grows 
according to 

i?2 _ SttG 

This dilutes the number of BHs at a rate 



{Prad + Pbh) ■ (24) 



( df{M,ty 



V dt ) --^f(^^4- (2^) 

\ / exp 

The total change in f{M, t) per unit time will result from the addition of the 5 contributions 
above. This fixes Pbh- 

= I iM M ^-M^ . (26) 

On the other hand, to obtain the change in T (or prad{T)) we impose energy conservation, 

diPradR") + dipBHR"") = '^Prad dR^ , (27) 
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where we have neglected the pressure of the BHs (they behave hke non-relativistic matter). 
The equation above imphes 

R R 

Prad = -4:Prad-^ " PBH - ^PBR-^ ■ (28) 

Notice that if the radiation and the BHs where decoupled, then one would have (\{pbhR^) = 
0, i.e., Prad = —4:Prad{R/ R)- This, however, is not the case since there is energy exchange 
between the two components. In particular, a variation in prad changes T, implying a change 
in the absorption/emission rate of the BHs and in Pbh- It is easy to see that 

d{pBHR^) =R^a dprad , (29) 

where we define a as 

a ^ ( ^] . (30) 



dprad 

Substituting this equation in (1271) we obtain 



R=cons. 



R^ (1 + a) dprad = Prad dR^ (31) 

The equation above expresses that BH evaporation can slow down the cooling of the 
radiation due to the expansion: as the universe expands, T drops, this leaves some BHs 
hotter than the plasma, so they evaporate and reheat the environment. Taking an interval 
where a is constant, eq. ( pTl) can be integrated to 

4 

Rq \ 1+" 



Prad ~ PradO J ' (32) 

If a ^ 1 BH evaporation would stop the change in the density and the temperature of the 
radiation due to the expansion. 



5 Black-hole dominated plasma 

We will now apply these equations to an initial configuration with only radiation (no BHs) 
at a given temperature T(0) = Tq. We find two generic scenarios depending on the value 
of Tq. Values closer to produce a larger number of BHs, that grow and absorb all the 
radiation before a Hubble time. Lower values of Tq imply a smaller number of BHs, that 
grow at basically constant temperature up to times of order H~^. In this section we will 
focus on the first case. We will describe the sequence of events using the toy model with 
n = 1, = 2, = 5, rric = 10 GeV, Md = 1 TeV and an initial temperature of Tq = 200 
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Figure 2: Distribution /(M,t) for U = 0.1, 1 GeV"^ and Tq = 200 GeV. Dashes describe the 
evolution neglecting BH-BH collisions. 

GeV. The critical BH mass (corresponding to Tbh = ^o) is Merit = 3.0 TeV. In Fig. 2 we 
plot the distribution f{M,t) for two values of t. We have included the distributions with 
and without the effect of collisions of two BHs to form a larger one. 

We can distinguish four phases in the evolution of the BH gas. 

1. In a first phase BH's of M > Merit are produced at constant rate (see Fig. 3) and 
absorb radiation of temperature T ^ Tq. As the number of BHs grows (see Fig. 2), 
collisions between two BHs to produce a larger one become important. This reduces 
the number of BHs and increses their average mass (in Fig. 3). At times around t ^ 2 
GeV~^ BHs start dominating the energy density and the temperature of the radiation 
drops (see Fig. 4). 

2. The drop in T stops the production of BHs in photon-photon collisions. In addition, 
the lighter BHs become hotter than the plasma, so they evaporate and feed prad- The 
evaporation reduces the number of BHs per unit volume, but the average BH mass 
grows: there is a continuous transfer of energy from the radiation and from the lighter 
(evaporating) BHs to the larger BHs of lower temperature. Once the BHs get a mass 
around 5 x 10'' GeV their radius stops growing. 

3. At t ~ 10^ GeV^^ the temperature of the radiation and of the BHs is similar, around 
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Figure 3: Number of BHs per unit volume and average BH mass as a function of time for 
To = 200 GeV. 

1.6 GeV. The (slow) mass growth of the heavier BHs is compensated by the decay of 
the lighter ones, with the radiation temperature basically constant. In this phase the 
energy density p = prad + Pbh is matter (BH) dominated. 

4. At times of order if ~^ ^ 2.6 x 10^^ GeV^^ the expansion cools the radiation. The BHs, 
of mass around 10^ GeV, decay fast (r ^ 4 x 10^ GeV~^) and the universe becomes 
radiation dominated. The lightest KK modes also decay fast (tkk ~ M^"^"'/mj?"'"" ~ 
10^ GeV~^), so only 4dim photons survive below T ^ 1 GeV. 

Two remarks are here in order. First, this is a toy model with only photons, gravitons, 
and BHs. In a more complete set up one should include baryons at T below 0.1 GeV (two 
types of matter, baryons and BHs, could coexist in models with lower values of mc/M^)). 
Second, this generic high T case, with Md > 5 TeV to avoid bounds from colliders and the 
inclusion of all the light standard model species (not just photons) at each temperature, 
could define a realistic set up. In these hybrid models it seems natural to obtain a plasma 
dominated by the standard particles at T ^ mc ^ ^qcd- The predictions for primordial 
nucleosynthesis would then be consistent with observations. 
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Figure 4: Temperature T of the radiation as a function of t. At larger values of t (up to a 
Hubble time ^ 10^^ GeV'^) T ^ 1.6 GeV. 



6 Radiation dominated plasma 

Let us now discuss the scenario with a low initial temperature and Pbh ^ Prad at any t. We 
take n = 1, = 2, gb = 5, = 10 GeV, = 1 TeV and Tq = 100 GeV. The critical BH 
mass (corresponding to Tbh = ^o) is Merit = 12 TeV, higher than in the previous case. As 
a consequence, the production rate of BHs colder than the plasma is much smaller. We can 
separate three phases in the evolution of this BH gas. 

1. At times below R-^ = 1.4 x 10^^ GeV"^ BHs are produced at constant rate in photon- 
photon collisions, with T ^ Tq. The number of BHs is so small that BH collisions can 
be neglected. All these BHs grow like the one discussed in Section 2. 

2. When the expansion cools the photons, BH production drops exponentially, whereas 
BH growth slows down. In Fig. 5 we plot the BH distribution f{M,t) after one and 
four Hubble times. The BHs reach a maximum mass around M = 10^^ GeV (there 
is just a 10% mass difference between 99% of the BHs) and Tbh ~ 1-6 GeV. This 
universe is always radiation dominated. At t ~ 10^^ GeV^^ the photon gas becomes 
colder than the BHs. 

3. In the last phase all the BHs evaporate in a time scale r ^ 10^^ GeV~^. 



13 



10-46 




M [GeV] 

Figure 5: Distribution /(M,t) for ti = H'^ = 1.4 x 10^^ GeV"^ and = 4//"^ with 
To = 100 GeV. 

This generic scenario is in principle consistent with primordial nucleosynthesis, since 
BHs are just a small fraction of the total energy density and do not change the expansion 
rate. In the particular case that we have discussed they decay when the plasma temperature 
is T ^ 0.01 GeV, but increasing the ratio Mu/rric one can obtain BHs that become 4 
dimensional and with a much longer lifetime. Their late decay could introduce distortions in 
the diffuse gamma ray background. In addition, in a more complete set up including baryons 
and structure formation they might work as seeds for macroscopic (primordial) BHs. 

7 Summary and discussion 

A transplanckian regime at accessible energies would have new and peculiar implications 
in collider physics and cosmology. What makes this regime special is that as the energy 
grows, softer physics dominate. For example, the production of a regular heavy particle at 
the LHC would provide an event with very energetic (hard) jets from its decay. In contrast, 
the production of a mini-BH would be seen as a high multiplicity event, with dozens of less 
energetic (soft) jets. Analogously, the production of massive particles or mini-BHs in the 
early universe would have very different consequences. The heavier the elementary particle, 
the shorter its lifetime, which tends to decouple these particles from low temperatures. For 
mini-BHs is just the opposite, heavier BHs are colder and live longer. In addition, if an 
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approximate symmetry makes a massive relic long lived, its late decay will produce very 
energetic particles. BHs would imply a much softer spectrum of secondaries [I9], with 
different cosmological consequences. 

In this paper we have analyzed the dynamics of a two-component gas with photons and 
BHs in an expanding universe. The system is characterized by the temperature T(t) of 
the radiation and the distribution f{M,t) of BHs. Our equations take into account BH 
production, absorption of photons, BH evaporation, and collisions of two BHs. We have 
discussed the two generic scenarios that may result from an initial temperature close to the 
fundamental scale of gravity M^. In the first scenario BHs empty of radiation the universe 
and dominate p before a Hubble time, whereas in the second case there are few BHs that 
grow at constant T up to t ^ H~^. 

We think that the work presented here is a necessary first step in the search for observable 
effects from these BHs. It could also lead us to a better understanding of the early universe 
at the highest temperatures. 
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